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I Abstract 

, We study the effects of heavy Majorana neutrinos on the Higgs sector 

■ of the MSSM via radiative corrections. We work within the SUSY con- 

text where the MSSM particle content is enlarged with right handed 
neutrinos and their corresponding SUSY partners, the sneutrinos, and 
where compatibility with neutrino data is required. We compute the 
one-loop corrections to the mass of the lightest MSSM CP-even neutral 
Higgs boson from Majorana neutrinos and their SUSY partners and 
assume a seesaw mechanism of type I for neutrino mass generation. A 
negative and sizeable Higgs mass correction of up to -5 GeV is found 
for a heavy Majorana mass of up to 10 15 GeV. This negative correc- 
tion can grow up to several tens of GeV if the soft SUSY breaking 
mass associated to their sneutrino partners is simmilarly heavy as the 
Majorana mass. 

Introduction 

The current experimental data on neutrino mass differences and neutrino 
mixing angles clearly indicate new physics beyond the so far successful Stan- 
dard Model of Particle Physics (SM). In particular, neutrino oscillations im- 
ply that at least two generations of neutrinos must be massive. Therefore, 
one needs to extend the SM to incorporate neutrino mass terms. 

We assume here the simplest version of a SUSY extension of the SM, the 
well known Minimal Supersymmetric Standard Model (MSSM), extended 
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by right-handed Major ana neutrinos and their SUSY partners, and where 
the seesaw mechanism of type I [1] is implemented to generate the small 
neutrino masses. For simplicity, we focus here in the one generation case 
which already provides interesting results and allows for a clearer illustration 
of the most relevant behaviour of the radiative corrections with the Majorana 
mass scale. The three generation case will be postponed for a later work. 

On the other hand, it is well known that heavy Majorana neutrinos, with 
rriM ~ 10 13 — 10 15 GeV, induce large LFV rates [2], due to their potentially 
large Yukawa couplings to the Higgs sector. For the same reason, radiative 
corrections to Higgs boson masses due to such heavy Majorana neutrinos 
could also be relevant. Consequently, our study has been focused on the 
radiative corrections to the lightest MSSM CP-even h boson mass, M/,, due 
to the one-loop contributions from the neutrino/sneutrino sector within the 
MSSM-seesaw framework. 

In the following we briefly review the main relevant aspects of the calcu- 
lation of the mass corrections and the numerical results. For further details 
we address the reader to the full version of our work [3], where also an 
extensive list with references to related works can be found. 



The basics of the MSSM-seesaw model 

The MSSM-seesaw model with one neutrino/sneutrino generation is de- 
scribed in terms of the well known MSSM superpotential plus the new rele- 
vant terms contained in: 



W = etj 



Y u H l 2 LPN - YiH\ V R 



+ ±Nm M N, (1) 



where mu is the Majorana mass and N = (u R , {vr) c ) is the additional 
superfield that contains the right-handed neutrino v R and its scalar partner 
vr. 

There are also new relevant terms in the soft SUSY breaking potential: 
Kjoft = + m\v* R i>R + {^ v A v H\v h v* R + m M B u i>RD R + h.c.) . (2) 

After electro-weak (EW) symmetry breaking, the charged lepton and 
Dirac neutrino masses can be written as 

mi = Yi vi , m D = Y u v 2 , (3) 

where V{ are the vacuum expectation values (VEVs) of the neutral Higgs 
scalars, with = v cos(sin)/3 and v = 174 GeV. 
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The 2x2 neutrino mass matrix is given in terms of mo and txim by: 



M" = ( ° mD ) . (4) 

Diagonalization of M v leads to two mass eigenstates, rij (i = 1, 2), which are 
Majorana fermions with the respective mass eigenvalues given by: 

m Vj n = ^ {mu T ^/ m 2 M + 4mf^ . (5) 

The mixing angle that defines the mass eigenstates is given by, 

m v m D 

tan 6* = = . (6) 

In the seesaw limit, i.e. when £ = <C 1 : 

rrij, = -m D £ + 0(m D £?) ~ — , m^v = ttim + 0(m D £) ~ m M , 

(7) 

and is small, therefore, i/ is made predominantly of z^l and its c-conjugate, 
(vl) C i whereas iV is made predominantly of vr and its c-conjugate, {vr) c - 

Regarding the sneutrino sector, the sneutrino mass matrices for the CP- 
even, M + , and the CP-odd, M_, subsectors are given respectively by 

M 2 - ( m | + m D + I M 1 cos 2Z 3 m D (A v - (j, cot f3±m M ) \ 
± I m D (A u - fi cot (3 ± mu) m 2 ^ + m 2 D + m 2 M ±2B v rriM J 

The diagonalization of these two matrices, M±, leads to four sneutrino mass 
eigenstates, hi (i = 1, 2, 3,4). In the seesaw limit, where rriM is much bigger 
than all the other mass scales the corresponding sneutrino masses are given 
by: 

m% +fi _ = m\ + cos 2/3 ^2m D (A v -n cot (3 - B v )£ , 

m % + ,N- = m M" ^ 2B v m M + m | + 2m D ■ ( 9 ) 

The mixing angles defining the sneutrinos mass eigenstates, 9±, are small in 
this limit and, therefore, z>+ and z>_ are made predominantly of and its 
c-conjugate, z>£, whereas N + and iV_ are made predominantly of vr and its 
c-conjugate, v* R . 



3 



Regarding the relevant interactions for the present work, there are pure 
gauge interactions between the left-handed neutrinos and the Z boson, those 
between the 'left-handed' sneutrinos and the Higgs bosons, and those be- 
tween the 'left-handed' sneutrinos and the Z bosons. All of these are com- 
mon to the MSSM. In addition, in this MSSM-seesaw scenario, there are 
interactions driven by the neutrino Yukawa couplings (or equivalently 
since Y v = (gm D )/(\/2M w sm(3)), as for instance gi WL v R = -fs^^fff > and 
new interactions due to the Majorana nature driven by myi which are not 
present in the case of Dirac fermions, as for instance g' hi>L y R = — %9 ^Mw'sinp ° • 

Finally, concerning the size of the new parameters that have been intro- 
duced in this model, in addition to those of the MSSM, i.e., tum, mD, 
A u and B u , there are no significant constraints. In the literature it is often 
assumed that tum has a very large value, mM ~ 0(1O 14-15 ) GeV, in order 
to get small physical neutrino masses \m u \ ~ 0.1 - 1 eV with large Yukawa 
couplings Y v ~ 0(1). This is an interesting possibility since it can lead to 
important phenomenological implications due to the large size of the radia- 
tive corrections driven by these large Yukawa couplings. We have explored, 
however, not only these extreme values but the full range for tum from the 
electroweak scale ~ 10 2 GeV up to ~ 10 15 GeV. 



The basics of the calculation 

In the Feynman diagrammatic (FD) approach the higher-order corrected 
CP-even Higgs boson masses in the MSSM, denoted here as and Mjj, 
are derived by finding the poles of the (h, i?)-propagator matrix, which is 
equivalent to solving the following equation [4]: 



P 



ml + T,hh(p 2 ) P 2 ~ m 2 H + £hh(p 2 



^hH(p 2 



0. 



(10) 



where rrihH are the tree level masses. The one loop renormalized self- 
energies, can be expressed in terms of the bare self-energies, 
S^(p 2 ), the field renormalization constants SZ^ and the mass counter 
terms Sm 2 ^, where (j) stands for h, H. For example, the lightest Higgs boson 
renormalized self energy reads: 



^hh{p 2 ) = ^hhip 2 ) + 5Z hh (p 2 - ml) 



5m 2 



hi 



[ID 



Regarding the renormalization prescription, we have used an on-shell 
renormalization scheme for Mz, My/ and Ma mass counterterms and T^, Th 
tadpole counterterms. On the other hand, we have used a modified DR 
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scheme (mDR) for the renormalization of the wave function and tan /3 . The 
mDR scheme is very similar to the well known DR scheme but instead of 
subtracting the usual A = f — 7b + log(4-7r) one subtracts A m = A — 
log(m\ I / /j^y^) , hence, avoiding large logarithms of the large scale uim- As 
studied in other works [5], this scheme minimizes higher order corrections 
when two very different scales are involved in a calculation of radiative 
corrections. 

Since we are interested in exploring the relevance of the new radiative 
corrections to the lightest Higgs mass from the neutrino/sneutrino sector, 
we will present here our results in terms of the mass difference with respect 
to the MSSM prediction. Consequently, we define, 

Am^:=Mf-M ft , (12) 

where M^ u denotes the pole for the light Higgs mass including the v/v cor- 
rections (i.e. in the MSSM-seesaw model), and Mh the corresponding pole in 
the MSSM, i.e without the v/v corrections. Thus, for a given set of input pa- 
rameters we first calculate Mh in the MSSM with the help of FeynHiggs [4], 
such that all relevant known higher-order corrections are included, and then 
we compute Am™ DR . The full one-loop neutrino/sneutrino corrections to 

the self-energies, Ei^, and entering in the evaluation of Am™ DR 
have been evaluated with the help of FeynArts [B] and FormCalc [7]. 



Analytical and Numerical Results 

We have fully explored the dependence with all the new parameters in [3] and 
found that the radiative corrections to the Higgs mass are mostly sensitive 
to niM, rn u , B u and nip,. We focus here on uim, rn^ and mp, leading to the 
largest corrections and address the reader to [3] for the complete study in 
terms of all the parameters. 

First, in order to compare systematically our predictions of the neu- 
trino/sneutrino sector in the MSSM-seesaw with those in the MSSM, we 
have split the full one-loop neutrino/sneutrino result into two parts: 

£ (p 2 ) | full = S (p 2 ) | g augc + S (p 2 ) | Yukawa , (13) 

where T>(p 2 )\ gauge means the contributions from pure gauge interactions and 
they are obtained by switching off the Yukawa interactions, i.e. by setting 
Y v = (or equivalently mo = 0). The remaining part is named here 
^(p 2 ) I Yukawa and refers to the contributions that are only present if Y v ^ 0. 
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In other words, this separation splits the full result into the common part 
with the MSSM, given by S(p 2 )| gaU g e , and the new contributions due to the 
presence of Major ana neutrinos with non vanishing Yukawa interactions, 
given by S(p 2 )| Yukawa- Thus, by comparing the size of these two parts, 
within the allowed parameter space region, we will localize the areas where 
Yukawa ^ ^ (p 2 ) I gauge, which will therefore indicate a significant depar- 
ture from the MSSM result. We see in Fig. [T]that the Yukawa part clearly 
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Figure 1: Left panel: Renormalized Higgs boson self-energy as a function 
of rriM- Here, p = 116 GeV. Right panel: Renormalized Higgs boson self- 
energy as a function of the external momentum p. In both panels, the two 
contributions from the gauge and Yukawa parts and the full result are shown 
separately. 

dominates for large tum ~ 10 14-15 GeV and that it shows a relevant external 
momentum dependence. Therefore, at large rriM, to keep just the Yukawa 
part is a good approximation, but to neglect the momentum dependence or 
to set the external momentum to zero are certainly not. In consequence, the 
alternative effective potential method will not provide a realistic result for 
the radiative corrections to the Higgs mass. 

Second, in order to understand in simple terms the analytical behavior 
of our full numerical results we have expanded the renormalized self-energies 
in powers of the seesaw parameter £ = raj^jmM'- 

s(P 2 ) = (fi(P 2 )) m „ + (S(P 2 )) mL , + (S(P 2 )) m , + .... (14) 

The zeroth order of this expansion corresponds to the gauge contribution 
and it does not depend on mp or mj,/. The rest of the expansion corresponds 
to the Yukawa contribution. The leading term of this Yukawa contribution 
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is the 0{m 2 D ) term, because it is the only one not suppressed by the Ma- 
jorana scale. In fact it goes as Y^M| W , where M§ w denotes generically 
the electroweak scales involved, concretely, p 2 , M§ and M\. In particular, 
the 0{p 2 m 2 D ) terms of the renormalized self-energy, which turn out to be 
among the most relevant leading contributions, separated into the neutrino 
and sneutrino contributions, are the following: 



y-mDR 



h- 



9 O 




VR 



-h + h---{ V--h 
- *- - 
VR 



m? D p 2 



g 2 p 2 m 2 D c 2 a 



+ 



g 2 p 2 m 2 D c 2 a 
M^M^s 2 ' 



(15) 



Notice that the above neutrino contributions come from the Yukawa inter- 
action ghv L v R , which is extremely suppressed in the Dirac case but can be 
large in the present Majorana case. On the other hand, the above sneutrino 
contributions come from the new couplings g'^ L ^ R , which are not present in 
the Dirac case. It is also interesting to remark that these terms, being ~ p 2 
are absent in both the effective potential and the RGE approaches. 

An useful formula for the Higgs mass correction has then be derived, by 
keeping just the dominant 0{m 2 D ) contribution to the Yukawa part of the 
renormalized self-energy: 



mDR 



2M h 
-g 2 m 2 D 



2M h 



\-2M 2 A cos 2 (a — P) cos 2 P + 2Ml cos 2 a 



128tt 2 M^M h sin 2 
-M| sin sin(a + /3) (2 (l + cos 2 0) cos a — sin 2/3 sin a)] . (16) 

We have checked that this simple expression reproduces very accurately the 



full numerical result, for the region of our interest with large mM~ 10 14 15 
GeV. 

With respect to the numerical results, figure [2] exemplifies the main fea- 
tures of the extra Higgs mass corrections Am™ DR due to neutrinos and 
sneutrino loops in terms of the two physical Majorana neutrino masses, 
and m v . For values of < 3 x 10 13 GeV and \m v \ < 0.1 — 0.3 eV the 
corrections to Mh are positive and smaller than 0.1 GeV. In this region, the 
gauge contribution dominates. In fact, the wider black contour line with 
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fixed Am™ = 0.09 coincides with the prediction for the case where just 
the gauge part in the self-energies has been included. This means that 'the 
distance' of any other contour-line with respect to this one represents the 
difference in the radiative corrections respect to the MSSM prediction. 

However, for larger values of mjv and/or \m u \ the Yukawa part domi- 
nates, and the radiative corrections become negative and larger in absolute 
value, up to values of -5 GeV in the right upper corner of Fig. These 
corrections grow in modulus proportionally to tum and m„, due to the fact 
that the seesaw mechanism imposes a relation between the three masses 
involved, m 2 D = \m u \mM- 




0.2 0.4 0.6 0.8 1.0 

\m v \ (eV) 

Figure 2: Contour- lines for the Higgs mass corrections from the neu- 
trino/sneutrino sector as a function of the physical Major ana neutrino 
masses, light \m v \ and heavy m^. The other parameters are fixed to: 
A v = B u = m- L = m R = 10 3 GeV, tan/3 = 5, M A = 200 GeV and 
^ = 200 GeV. 

Finally, we plot in Fig. [3l the contour-lines for fixed Am" DR in the 
less conservative case where mp, is close to tum- These are displayed as a 
function of \m v \ and the ratio m^jmu- m-M is fixed here to the reference 
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value, m,M = 10 GeV. For the interval studied here, we see again that 
the radiative corrections can be negative and as large as tens of GeV in 
the upper right corner of the plot. For instance, Am™ DR = —30 GeV for 
m,M = 10 14 GeV, \m u \ = 0.6 eV and m^/ruM = 0.7. 




0.2 0.4 0.6 0.8 1.0 



\m v \ (eV) 

Figure 3: Contour-lines for the Higgs mass corrections from the neu- 
trino/sneutrino sector as a function of the ratio m^/mM and the light- 
est Majorana neutrino mass \m u \. The other parameters are fixed to: 
m M = 10 14 GeV, A v = B u = m L = 10 3 GeV, tan/3 = 5, M A = 200 GeV 
and fj, = 200 GeV 



Conclusions 

We have used the Feynman diagrammatic approach for the calculation of the 
radiative corrections to the lightest Higgs boson mass of the MSSM-seesaw. 
This method does not neglect the external momentum of the incoming and 
outgoing particles as it happens in the effective potential approach. We 
have performed a full calculation, obtaining not only the leading logarithmic 
terms as it would be the case in a RGE computation but also the finite 
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terms, that we have seen that can be sizable for heavy Majorana neutrinos 
(10 13 — 10 15 GeV) and the lightest neutrino mass within a range inspired by 
data (0.1 — 1 eV). For some regions of the MSSM-seesaw parameter space, 
the corrections to are substantially larger than the anticipated LHC 
precision (~ 200 MeV) [8]. Specifically they can be negative and up to -5 
GeV if rriM ~ 10 15 GeV and \m u \ ~ 1 eV and up to minus tens of GeV if in 
addition is also large and of similar size to m^j. 
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